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Abstract 

We consider a 2D isotropic Fermi liquid with attraction in both s and d 
channels and examine the possibility of a superconducting state with mixed s 
and d symmetry of the gap function. We show that both in the weak coupling 
limit and at strong coupling, a mixed s + id symmetry state is realized in a 
certain range of interaction. Phase transitions between the mixed and the 
pure symmetry states are second order. We also show that there is no stable 
mixed s + d symmetry state at any coupling. 
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I. INTRODUCTION 



The question of the order parameter symmetry is one of the central issues of high- 
temperature superconductivity. There is a general consensus that the superconducting gap 
is highly anisotropic, but whether the gap has a particular symmetry under rotations is still 
a matter of debate. A number of experiments on YBCO-123 are roughly consistent with 
the d— wave symmetry JT] for which the most natural source is the exchange of magnetic 
fluctuations 0, but some experiments, e.g., photoemission studies on BSCCO-2212 |3|,|J as 
well as c-axis Josephson tunneling experiments on YBCO |Q, are inconsistent with the pure 
d-wave, but more consistent with ans + cf state. In principle, the presence of orthorhombic 
distortion in, e. g. , BSCCO-2212 || in itself guarantees that an otherwise o?-wave supercon- 
ducting gap will have an admixture of the s-wave [0. However, the superconducting state 
may be a mixture of s and d components even in the absence of an orthorhombic distortion. 

A superconducting state with a mixed s + d symmetry of the gap was first discussed in 
||, and the s + id state in 0. The mixed symmetry state at intermediate doping levels 
was also found in variational Monte Carlo studies of the t — J model ||10|| . An alternative 



possibility of symmetry mixing due to interplane coupling was proposed in Ref. |TT|. More 
exotic mixed symmetry states have also been suggested . 

Very recent work has shown [nj that the extent of s-wave admixture is a strong function 
of t', the second neighbor hopping, which varies a good deal from one high-T c material to 
another. This strongly suggests that the question of s-d mixing should be looked at as a 
function of hole doping, and that this must be done in each high-T c material separately. 

The variational Monte Carlo calculations indicate that at the doping levels which favor 
mixed symmetry states, the ground state energy is roughly independent of the relative 
phase 9 of s- and d-wave components. This relative phase is of great importance, since 
only 6 = and 9 = n states can have gap nodes. The experiments on combinations of 



Josephson junctions on YBCO-123 [15| appear to rule out a relative phase of ir/2 if the 
s-wave amplitude is more than about 10% of the d-wave. Recent photoemission work as 
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a function of hole doping [[|, ]T6| indicates that 6 = and that the relative amplitude of 



s-wave and d-wave depends on temperature and hole doping level. 

In view of both the experimental situation, which appears to indicate that s-d mixing 
is possible in some systems and the theoretical situation where the relative phase is not 
well-determined, additional understanding of the physics of this phase is needed. In this 
paper, we address the issue whether one can obtain the s + d mixed state in calculations on 
a simple but fairly general model. The answer we obtain is negative - we found that for any 
coupling, the only possible mixed symmetry in this model is s + id. 

We consider a model of an isotropic two dimensional Fermi liquid with attractive interac- 
tion in both s and d channels. We assume that both interactions are frequency-independent 
in a frequency range bounded by the cut-off frequency u> c , and are zero for \u\ > uj c . Obvi- 
ously, when only one of the two interaction channels is present, the ground state is described 
by the corresponding pure symmetry gap function. When both interaction channels are 
present, their competition will lead to either one of the two pure symmetry superconducting 
states, or a mixed state, where the gap function contains both the s and the d harmonics. 

In the next section we will consider the weak coupling case, where one can use the BCS 
formalism. We will show that the transition between s and d symmetries occurs via an 
intermediate phase with mixed s + id symmetry. The two phase transitions between the 
pure and the mixed states are second order. In Section III we consider the case of arbitrary 
coupling in the framework of the Eliashberg theory. We will show that there always exists 
a range of relative strengths of the s and d interactions where an s + id solution exists. The 
analysis of the s + d mixed state is more complicated. However, we can show that at least 
in both the weak and the strong coupling limits the s + d mixed state does not occur. Our 
conclusions are summarized in Section IV. As an aside, in the Appendix we also present 
few simple results for the thermodynamics of a d-wave superconductor in the weak coupling 
limit, which to the best of our knowledge have not been published anywhere else. The main 
feature is that the ratio of the superconducting gap to the transition temperature, 2A/T C , 
for d-wave is 4.28, larger than 3.06 for the s-wave . 



II. WEAK COUPLING 



In this section we will consider the case when the coupling is weak in both interaction 
channels. In this case the BCS theory is valid, and the gap equation assumes the following 



form [17 1 
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V m = V a (k, k') + V d (k, k') cos(2(0 - 0')) 
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is the interaction, which contains both s and d harmonics. The angle is defined in our two- 
dimensional model as = tan" 1 (k y /k x ). Consider now a trial mixed state with arbitrary 
phase difference between the s-wave and the d-wave components of the order parameter: 



AO) = A s + e i9 A d cos: 



(4) 



Separating the real and imaginary parts of this equation, and also s and d components, we 
obtain a set of three independent equations 
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where V s = V s (kp,k F ) and V d = Vd(kp, hp)- It is straightforward to see that if both A s 
and A d are finite, the set ([5] - 0) can be simultaneously satisfied in only two cases 9 = or 



9 = 7r/2, leading to s + d or s + id respectively. Thus the weak-coupling theory gives the 
same restricted set of possibilities for the internal phase angle that Ginzburg-Landau theory 
offers flEfl . Below we consider these two cases separately. 



A. s+d state 



In this case 9 = 0, and equations (|5|) and (^) become : 

o 27r 4^/e 2 + A 2 + 2A s A d cos 20 + A 2 d cos 2 20 

Jo 2tt J-u, e 4 J e 2 + A 2 + 2A s A d cos 20 + A 2 cos 2 20 



Performing the integration over frequency and doing standard manipulations, we obtain 

(l-^i)a = g s f(a) (11) 
9d 

where g s = — K,iV(0)/4, g d = — V d N{0)/A, a = A s /A d , and the function /(a) is given by 

f(a)= / — (2a cosx — l)(cosx + a) log(a + cosx) 2 . (12) 
jo 2n 

The graphical solution of Eq.([0]) is shown in Fig. [1]. It is easy to see that in the limit of 
a — > oo, /(a) ~ a, while in the limit of a — > 0, /(a) ~ —a/2. If we start out with g s = 0, 
then, naturally, the only solution is a = 0, i.e. pure c?-wave. As g s increases, the slope of 
the straight line on Fig. [TJ decreases, and at 

9? = • (13) 

the lines first cross at a = oo. If we increase g s even further, we find three solutions: at 
a = 0, a = oo, and at some finite a m , which corresponds to a mixed state. As g s increases 
further, a m decreases and becomes zero at 

( 2) = __9d (14) 
2-g d /2 1 ; 

For larger g s , there exists only one solution: a = oo which corresponds to a pure s-wave 

state. 



We see therefore that d— wave solution exists at < g s < while the s— wave solution 
exists at g s > g^ 1 '. The key point is that g^ > g^\ such that there is an intermediate region 
g^ < g s < g^p where both pure solutions exist together with the s + d solution. To verify 
which solutions are stable, we computed the second derivatives of the energy and found that 
the two pure solutions are stable in the intermediate region while the s + d state actually 
corresponds to a maximum rather than a minimum of energy. Clearly then, s + d state is 
unstable; if it was the only mixed state allowed, then the system would simply undergo a 
first-order transition between the two pure states with a region where hysteresis is possible 
between gjp and gf\ 



B. s+id state 

In the case of 6 = ir/2, we follow the same procedure. Now the coupled gap equations 
have the following form 

Ad = f d JL f A rf cos 2 20' 

Jo 2ir J- a , J e 2 + A 2 + A 2 cos 2 20 



A, = t f r dig , A - . (16) 

Jo 2vr J-u c e 2 + A 2 + A 2 cos 2 2( h 



Integrating over frequency and doing standard manipulations we obtain that the mixed state 
exists if 

1 _ 2 Jl = 9s ( I + a 2 - aVTT^) . (17) 

As before, a = A s /A d . 

As we start out with g s = 0, Eq. ( p!7|) has no solution and the gap has pure d symmetry 
(a = 0). However, contrary to the previous case, here a solution of fll7|) first appears at 

flC") = i* (18) 

for the same a = 0. As we increase g s from g^ ll \ a and therefore A s increases continuously, 
satisfying 



g s /2 + 2g s /g d -l 
a = , (19) 

pg s {l-2g s /g d ) 

and becomes infinite at g^ = g&j2. Clearly, in this situation, we have a second order phase 
transition from pure d to a mixed s + id symmetry state at g s = g^ % \ and a second order 
transition from a mixed state to pure s state at g s = g[ 2 ^ > g^. In other words, the 
two pure states, which are stable with respect to s + d mixture, are in fact unstable (for 
corresponding g s ) with respect to s + id mixture, and in between g^ and gf % \ the s + id 
state is the equilibrium state of a system (see Fig.^J) . 

III. STRONG COUPLING 

In order to be certain that our results are not an artifact of the weak-coupling approxi- 
mation, we perform the calculations in the strong coupling regime. We follow the Eliashberg 
formalism JTB||TP| ] at zero temperature. We assume that the frequency cutoff uj c <C ep, so 
that vertex corrections can be neglected according to Migdal theorem |^D| . 



In the Eliashberg approach, one preserves the frequency dependence of the gap, and 
substitutes the full quasiparticle Green function in the gap equation. In explicit form, the 
equations are 



l? Q 2 {k', u') + A 2 {k', uj>) + £ 2 (k>) 



Q(k, uj)=uj- / Y V tl :, ^ )r- L =r 1 ^ , 20 

V ^ J 2lx 9 n 2 (k' 1 uj') + A 2 (k',uj l )+e(k') 

where £(<f) is the renormalized single particle energy, £(g) = Vp(q — Pf), A(q,u) is the 
anomalous part of the self-energy, and Q(q,Lj) is the antisymmetric in uo part of the normal 
self-energy. As uj c <C e^, the integration over \k'\ is confined to a region near the Fermi 
surface and can be substituted by the integration over £, as in BCS theory. This integration 
is straightforward, and yields 

, /**W2 , d(p (g s + g d cos2(p)A(u', 
Aw = / duo 



u,-u /2 JO 27T ^1^(^0)12+ | A(U',4>)\< 



Jo ^ ^J\n(uj',(f))\' 2 + \A(u',(l))\ 2 
We now again consider the two mixed states separately. 



A. s + id state 

Consider first a mixed s + id state. For this symmetry, the angular decomposition of the 
self-energy functions yields 

A(w, 0) = A s (u) + iA d (u) cos(20) (22) 
0) = tt s (uj) + itt d (uj) cos(20) . (23) 

Accordingly, Eqs.(^) can be decomposed into four equations 

a„m = r- n ^ r a ± , (24) 

Ju-u»/2 Jo 2?r ^Ja 2 s (ui') + n 2 d (ui') + A 2 (w') + A 2 .^') COS 2 2<* 
i.M = P"' 2 ^' f ^ , 9 - A ' ( "'> (25) 

Ju-ua/2 JO 27T y^^/) + Q2^/) + £2^/) + A 2 (cj') COS 2 20 

^ = p o/2 *s r - 1 sm^^ 

Ju-u /2 Jo 2tt ^(u/) + n 2 d (oo') + A 2 (u/) + A 2 {uo') cos 2 20 
= W + ^ ; ££±2 (27) 

A,-u, /2 ^0 27T ^Q^/) + fljj^/) + £2( w /) + A 2 (cj') COS 2 20 

Equation (^) is homogeneous in For weak coupling its only solution was Q d = 0. In 

principle, at strong coupling, there is a chance that above some threshold, there exists a 

solution with a non-zero f^. We will not consider this rather exotic possibility and assume 

that the stable solution of Eq.(p6|) corresponds to = for all couplings. 

We now follow the same approach as at weak coupling and look for the transition points 

between pure s and mixed s + id, and pure d and mixed s + id states. In the former case, 

we linearize the above set of equations around A^ = and obtain 

. r"+"o/2 [2* d<j) g d A d (u') cos 2 20 
AAuj) = did 



u-wo/2 Jo 2tt ^n 2 s (LU>) + A 2 s (lu') 

~ , , ru+ua/2 r2n d<h q s A s (u') , , 

A s (u) = dJ\ — . ys sK ' (28) 

Jv-uo/2 JO 27T J&ftf) 



Eqs.([28|) are obviously satisfied when = gd/2, the same as for weak coupling. 

Now consider the transition from the (i-wave state into the mixed state. Linearizing 
Eqs.fl2|)-(p7D with respect to A s we get 



Jo 2tt 2 . kWuA + te(u/) cos 2 



+ r +UJ ° /2 duj' F* — &Az(^)cos40 ^ 
J"-"o/2 Jo 2tt 2y/£l 2 8 (u') + A 2 (c/) cos 2 2cf) 

A, M = dJ r d ± , tht2 (30) 

Jv-uo/2 Jo 2n / Q 2 (u')+h 2 d (u>) cos 2 20 



We first observe that if the second term in Eq. ( |2"9"D were absent, the d— wave solution 
would become unstable at gf^ = gd/2 = gf^- We now show that this second term yields 
a negative correction to the first term independent of what A^(a/) and Q s (u') are. Indeed, 
let's perform angular integration first. For the first term in (p9"D, the integrand (apart from 
gdAd(u')) is positive, and the integration yields a positive result. For the second term, we 
have to evaluate 



'o y/a + b cos 
where a > b > 0. Doing simple manipulations, we obtain 

rn / 2 , , / COS0 COS 



Jo 



(32) 



\ \J a + b cos y/a — b cos / 
Since cos0 is positive when < < tt/2, then I is negative. This simple argument shows 
that the second term in (p9|) effectively reduces the value of ^ to g e d < g^. Clearly then, 
the critical value g^ = g d ' j2 for the transition between the pure d state and the s + id 
mixed state will be smaller than gd/2 = g( 2t \ 

The above consideration shows that for arbitrary interaction strength, exists a range 
°f gs/gd where neither of the pure states is stable. We further expanded near each of the 
transition points up to qubic terms in either A s or Ad, and indeed found a nonzero solution 
for this intermediate range of couplings which implies that in the intermediate region, the 
gap has an s + id symmetry. It turns out therefore that the phase diagram at weak and 



strong couplings is essentially the same. It is nevertheless interesting that the strong coupling 
corrections tend make the mixing of s-wave into a predominantly d-wave more favorable. 



B. s + d state 

Now we study whether it is possible to obtain the s + d mixed symmetry state in equilib- 
rium. This case is significantly more complicated, because unlike the s + id case, the square 
of the gap function now contains a term which is linear in both s and d components. 

However, we will show that at least in the limit of very strong coupling, there is no region 
of s + d symmetry. Indeed, consider the transition between the pure s and the mixed s + d 
states. Eliashberg equations linearized in A^ and Qd read 

Ju-m/2 Jo 2jr y'fiZfu') + A 2 d (ui') 

n.M = « + r°< 2 ^ r (3 4) 



L-uo/2 Jo 2ir y/Wftc/) + A 2 (a/) cos 2 
g d A d {u'){A d {u')A 8 {u)') + Q d (u')Q s (uj')) cos 2 
2(f) 2 (u/) + A 2 V)) 3/2 



(35) 



JO 27T ^2(^) + A^COS 2 

g d A d {u'){A d {uj')A s {u') + cos 2 



(36) 



2(0 2 (a,') + A 2 (^)) 3/2 

We now show that for g s ^> 1 and < u , A s (lj) ^> fi s (c<j). Indeed, suppose that this is 
true. Then it is easy to see that in the region of frequencies we are interested in, A s (u) is 
frequency independent and equal to g s u>o, while fl s (u) is linear in uj. Letting VL s {uj) = \u, 



we can rewrite the inequality as g s ^> A. Then, equation fl34]) becomes 



q s X JdJ . 

A = 1 + — / . . (37) 

u J^ /2 ^AV 2 + A 2 



Solving ( P?[ ) we obtain for g s ^> 1 
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A = 2g 2 J 3 « g s 



(38) 



thus justifying the assumption that A (a;) 3> £l s (u). Furthermore, in this limit, Eq.(p5[) 



reduces to 




-0/2 \ 2 uJ 2 duJ 



(39) 



0,0/2 (A 2 , + A 2 a; 2 ) 3 /2 



Using (|38 ) we finally obtain that 



(40) 



We see therefore that in the strong coupling limit, the hypothetical transition between a 



couplings, both the pure s-wave and the mixed state clearly must be unstable with respect 
to the pure (i-wave state. A similar analysis shows that the transition between the pure s 
and the mixed state occurs at g s = g^ 1 ' ~ g d . As a result, we again have g^ > gf \ which, 
just as in the weak coupling case, implies that there is no region of mixed s + d symmetry. 



We have studied in this paper a two-dimensional isotropic Fermi liquid with attractive 
interaction in both s and d channels. We considered the weak coupling limit and also applied 
the Eliashberg formalism at strong coupling. 

The phase diagram of the superconductor turns out to be the same at weak and strong 
coupling. It displays a region with a mixed s + id symmetry gap when the coupling strengths 
in the two channels are of the same order of magnitude. The phase transitions between the 
mixed state and the pure s and d states are second order. 

We have shown that in the weak and in the strong coupling limits a mixed s + d state 
does not occur. Intuitively, this can be interpreted as the propensity of the system to choose 
the state in which the amplitude of the gap function has the largest value. This is also 



mixed s + d and a pure s states occurs at g s = g^ gd- However, for this ratio of 



IV. SUMMARY 
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in agreement with the Ginsburg-Landau considerations [|Tl],|nj, which suggest that in the 



absence of orthorhombic distortion, the s + id state has lower energy than the s + d state. 

Indeed, the model we considered is oversimplified: the two-dimensional isotropic Fermi 
liquid captures some of features of the high-Tc materials, however, the isotropic system is 
never close to a magnetic or a metal-insulator transition. Our analysis of the s + d versus 
s + id question therefore does not include such effects. In this sense, what we have shown 
here is that if the s + d state is realized in real materials, some nontrivial physical effect of 
the proximity to these phase transitions is likely to be the cause. 
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VI. APPENDIX 

In this Appendix, we present weak coupling calculations for the thermodynamics of a 
pure ci-wave superconductor. 

The gap equation at zero temperature is 

rwo/2 f2n ^0 COS 2 26 , . 

1 = 9d / de — . . (41) 

J-uo/2 Jo 27T ^ e 2 + A 2 cog 2 20 

After simple manipulations we arrive at 



A, = ^expf-i). (42) 
Ve V 9dJ 



At the same time, the transition temperature is given by [17] 



n \ 9d, 
where log 7 = C ~ 0.577 is Euler's constant. Then 



r c = 2^exp(-I) , (43) 
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2A r , 



47T 



4.28 . 



Note, this ratio for s-wave is 3.06. 

Finally, starting with the gap equation at finite temperature 



(44) 



9d 



w /2 
-wo/2 



de 



^ d(j) cos 2 2<j> 



o 2vr e 2 + A 



(45) 



and performing standard expansions at T <C T c and T C — T<^T C []T7 we find the following 
expressions for the temperature dependence of the superconducting gap 



A d (T) = A d (0) 



ji\ 3 

1-0.371- 



T <T C 



A d (T) = 1.65A d (0Wl - T/T c T c -T < T c 



(46) 
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FIGURES 

FIG. 1. Graphical solution of Eq. (12) is given by the intersection of the curve f(a) and the 
straight line. 

FIG. 2. The location of the phase boundaries for the trial s + d state. The critical point of the 
transition from s + d to s occurs inside the d-phase and vice versa, meaning that there is no region 
of s + d mixed phase. 

FIG. 3. The phase diagram of the superconductor at zero temperature depending on the ratio 
of coupling strengths in the two channels. 
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